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Abstract 
l/^ ' 

^ ' Filtration combustion is described by Laplacian growth without surface ten- 

, sion. These equations have elegant analytical solutions that replace the complex 

' integro-differential motion equations by simple differential equations of pole mo- 

QQ ' tion in a complex plane. The main problem with such a solution is the existence 

, of finite time singularities. To prevent such singularities, nonzero surface ten- 

sion is usually used. However, nonzero surface tension does not exist in filtration 
] combustion, and this destroys the analytical solutions. However, a more elegant 

approach exists for solving the problem. First, we can introduce a small amount of 
pole noise to the system. Second, for regularisation of the problem, we throw out 
all new poles that can produce a finite time singularity. It can be strictly proved 
that the asymptotic solution for such a system is a single finger. Moreover, the 
' qualitative consideration demonstrates that a finger with i of the channel width is 

, statistically stable. Therefore, all properties of such a solution are exactly the same 

as those of the solution with nonzero surface tension under numerical noise. The 
solution of the ST problem without surface tension is similar to the solution for the 
equation of cellular flames in the case of the combustion of gas mixtures. 

Keywords. Saffman - Taylor Problem, final time singularity, Laplacian Growth, Hele- 
Shaw cell, zero surface tension. Filtration Combustion, pole solution 

1 Introduction 

The problem of pattern formation is one of the most rapidly developing branches of 
nonlinear science today llt ilSJ . 
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The 2D Laplacian growth equation describes a wide range of physical problems, 
for example, filtration combustion in a porous medium, displacement of a cold liquid 
in a Hele-Shaw channel by the same liquid that is heated or of a hot gas in a Hele- 
Shaw channel by the same gas that is cooled, and solidification of a solid penetrating 
a liquid in a channel \2u 4i- This equation has an elegant analytical solution in the 
form of poles. However, such an equation can lead to the appearance of final time 
singularities. To prevent these singularities and to regularise the problem, a term con- 
taining the surface tension is usually introduced into the equation describing Laplacian 
growth. Unfortunately, in the presence of such a surface tension term, obtaining an 
analytical solution in the form of poles becomes impossible. In addition, it is usually 
assumed that the surface tension explains the occurrence of an asymptotic solution in 
the form of a finger with half of the channel width. This asymptotic behaviour is also 
observed in experiments. In this paper, the mathematical mechanism of the regular- 
isation is introduced. It makes it possible to avoid final time singularities, results in 
desirable asymptotic behaviour in the form of a finger with half of the channel width, 
and maintains the analytical solution in the form of poles. Maintenance of the analyt- 
ical character of the solution is very important - it makes it possible to easily analyse 
2D Laplacian growth solutions and to qualitatively or quantitatively explain the be- 
haviour The author sincerely hopes that this paper will play the same role for the 2D 
Laplacian growth equation as the paper LI 7, 1 did for the theory of gaseous combustion 
of pre-mixed flames. The ifTTl analytical solution and its asymptotic behaviour have 
given a push to development of the theory of gaseous combustion of pre-mixed flames 
and have made it possible to qualitatively or quantitatively explain the behaviour of a 
front of pre-mixed flames lfT2HT6l[T8l . 

Prof. Matkowsky, in collaboration with Dr Aldushin, considered planar, uniformly 
propagating combustion waves driven by the filtration of gas containing an oxidiser, 
which reacts with the combustible porous medium through which it moves ||2}S1. 
These waves were typically found to be unstable with respect to hydrodynamic per- 
turbations for both forward (coflow) and reverse (counterflow) filtration combustion 
(FC), in which the direction of gas flow is the same as or opposite to the direction of 
propagation of the combustion wave, respectively. 

The basic mechanism leading to instability is the reduction of the resistance to flow 
in the region of the combustion products due to an increase of the porosity in that re- 
gion. Another destabilising effect in forward FC is the production of gaseous products 
in the reaction. In reverse FC, this effect is stabilising. In the case in which the planar 
front is unstable, an alternative mode of propagation in the form of a finger propagat- 
ing with constant velocity was proposed. The finger region occupied by the combustion 
products is separated from the unburned region by a front in which chemical reactions 
and heat and mass transport occur 

In the papers [2 -4] of Prof. Matkowsky and Dr Aldushin, it was shown that 
the finger solution of the combustion problem can be characterised as a solution of 
a Saffman-Taylor (ST) problem originally formulated to describe the displacement of 
one fluid by another having a smaller viscosity in a porous medium or in a Hele-Shaw 
configuration. The ST problem is known to possess a family of finger solutions, with 
each member characterised by its own velocity and each occupying a different fraction 
of the porous channel through which it propagates. The scalar field governing the 
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evolution of the interface is a harmonic function. It is natural, then, to call the whole 
process Laplacian growth. 

The mathematical problem of Laplacian growth without surface tension exhibits a 
family of exact analytical solutions in terms of logarithmic poles in the complex plane. 

The main problem with such a solution is existing finite time singularities. To 
prevent such singularities, nonzero surface tension usually is used ( |5 -9|). The surface 
tension also results in a well-defined asymptotic solution: only one finger with half of 
the channel width. In addition, the other terms can be used for regularisation (see 
mOlll II and references therein). 

The solution of the ST problem without surface tension is similar to the solution 
for the equation describing cellular flames in the case of combustion of gas mixtures 
ifTzHTSll . Indeed, in both cases, solutions can be transformed to the set of ordinary 
differential equations. This set describes the motion of poles in the complex plane. 
Applying nonzero surface tension to the ST problem destroys this elegant analytical 
solution. 

It must be mentioned that the filtration combustion and the gaseous combustion 
in pre-mixed flames are features of different physics; the equation of 2D Laplacian 
growth and the equation describing the Mihelson-Sivashinsky feature use completely 
different mathematics. Moreover, whereas the equation for the Mihelson-Sivashinsky 
poles involves trigonometric functions, the equation for the 2D Laplacian growth poles 
involves logarithmic functions. The analogy here is not "half-baked" but rather is deep. 
Indeed, the very different and complex integro-differential equations have a simple 
analytical solution in the form of poles. Moreover, even the behaviour of these poles is 
similar. 

Another problem is the fact that surface tension may not be introduced for the math- 
ematical problem considered by Saffman and Taylor involving filtration combustion in 
a porous medium ||2}S1. Here, the zone of chemical reaction and diffusion of heat and 
mass shrinks to an interface separating the burned region from the unburned region. In 
all these problems, there is no pressure jump at the interface, so surface tension may 
not be introduced. Thus, the Saffman-Taylor model arises not only as the limiting case 
of zero surface tension in a problem in which surface tension enters the problem in 
a natural way but also in other situations in which the introduction of surface tension 
makes no sense. Another such problem is that of the solidification of a solid penetrat- 
ing a liquid in a channel. It is reasonable to expect that the selection may be affected 
by introducing a perturbation other than surface tension that is relevant to the specific 
problem under consideration. For example, in the combustion problem, the effect of 
diffusion as a perturbation might be considered. Here, the effect of diffusion is similar 
to that of surface tension in the fluid displacement problem |I2}S1. 

Therefore, we need to look for a solution without introducing surface tension using 
different methods for regularisation. 

Criteria were proposed ( |i3j and [4]) to select the correct member of the family of 
solutions (one finger with half of the channel width) based on a consideration of the ST 
problem itself, rather than on modifications of the problem. A modification is obtained 
by adding surface tension to the model and then taking the limit of the vanishing surface 
tension. 

Unfortunately, it is not clear from the papers why Laplacian growth without surface 
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tension gives this asymptotic solution (one finger with half of the channel width) that 
satisfies the identified criteria. 

However, a more direct way exists to solve the problem. First of all, we can intro- 
duce a small amount of noise to the system. (The noise can be considered a pole flux 
from infinity.) Second, for regularisation of the problem, we throw out all new poles 
that can produce a finite time singularity. It can be strictly proved that the asymptotic 
solution for such a system is a single finger Moreover, the qualitative consideration 
demonstrates that a finger with ^ of the channel width is statistically stable. Therefore, 
all properties of such a solution are exactly the same as for the solution with a nonzero 
surface tension under numerical noise. 

The rest of the paper is organised as follows. We begin by presenting arguments 
about Saffman-Taylor "finger" formation with half of the channel size (Section |2]l. 
Next, Section[3]describes asymptotic single Saffman-Taylor "finger" formation without 
surface tension. Finally (Section|4]i, we provide a summary and conclusions. 

2 Saffman-Taylor "finger" formation witli lialf of tlie 
channel size 

The case of Laplacian growth in the channel without surface tension was considered in 
detail by Mark Mineev-Weinstein and Dawson |19|. In this case, the problem has an 
elegant analytical solution. Moreover, they assumed that all major effects in the case 
with vanishingly small surface tension may also occur without surface tension. This 
would make it possible to apply the powerful analytical methods developed for the no 
surface tension case to the vanishingly small surface tension case . However, without 
additional assumptions, this hypothesis may not be accepted. 

The first objection is related to finite time singularities for some initial conditions. 
Actually, for overcoming this difficulty, a regular item with surface tension was intro- 
duced. This surface tension item results in loss of the analytical solution. However, 
regularisation may be carried out much more simply - simply by rejecting the initial 
conditions that result in these singularities. The second objection is given in work by 
Siegel and Tanveer p20l. There, it is shown that in numerical simulations in a case with 
any (even vanishingly small) surface tension, any initial thickness "finger" extends up 
to i the width of the channel. The analytical solution in a case without surface ten- 
sion results in a constant thickness of the "finger" equal to its initial size, which may 
be arbitrary. Siegel and Tanveer, however, did not take into account the simple fact 
that numerical noise introduces small perturbation to the initial condition or even dur- 
ing "finger" growth, which is equivalent to the remote poles, and with respect to this 
perturbation, the analytical solution with a constant "finger" 

It was shown by Mark Mineev-Weinstein 1211 that similar pole perturbations for 
some initial conditions, can be extended to the Siegel and Tanveer solutions. This pos- 
itive aspect of the paper 121] was mentioned by Sarkissian and Levine in their Com- 
ment ll22l . In summary, it is possible to determine that to identify the results with and 
without surface tension, it is necessary to introduce a permanent source of the new re- 
mote poles; the source may be either external noise or an infinite number of poles in 
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an initial condition. Which of these methods is preferred is still an open question. 

In the case of flame front propagation, it was shown [12-16| that external noise is 
necessary for an explanation of the flame front velocity increase with the size of the 
system: using an infinite number of poles in an initial condition cannot give this result. 
It is interesting to know what the situation is in the channel Laplacian growth. One 
of the main results of Laplacian growth in the channel with a small surface tension 
is Saffman-Taylor "finger" formation with a thickness equal to ^ the thickness of the 
channel. To use the analytical result obtained for zero surface tension, it is necessary 
to prove that formation of the "finger" also takes place without surface tension. 

In our teamwork with Mark Mineev-Weinstein ||23l , it was shown that for a finite 
number of poles at almost all allowed (in the sense of not approaching finite time sin- 
gularities) initial conditions, except for a small number of degenerate initial conditions, 
there is an asymptotic solution involving a "finger" with any possible thickness. Note 
that the solutions and asymptotic behaviour found in i23i for a finite number of poles 
are an idealisation but have a real sense for any finite intervals of time between the ap- 
pearance of the new poles introduced into the system by external noise or connected to 
an entrance to the system of remote poles of an initial condition, including an infinite 
number of such poles. The theorem proved in ||231 may again be applied for this final 
set of new and old poles and again yields asymptotic behaviour in the form of a "fin- 
ger", but the thickness is different. Thus, introduction of a source of new poles results 
only in possible drift of the thickness of the final "finger" but does not change the type 
of solution. 

It should be mentioned that instead of periodic boundary conditions, much more 
realistic "no flux" boundary conditions may be introduced ll24l . (This paper repeats 
the result for single finger asymptotic behaviour already proved in the papers ||231 ). 
This result forbids a stream through a wall, which inserts additional, probably useful 
restrictions on the positions, number, and parameters of new and old poles (explaining, 
for example, why the sum of all complex parameters ai for poles gives the real value 
a for the pole solution (5) in [21 1). However, this does not have an influence on the 
correctness and applicability of the results and methods proved in ||23l . 

In EJK Mark Mineev-Weinstein gives "proof that steady asymptotic behaviour 
for Laplacian growth in a channel with zero surface tension is a single "finger" with 
a thickness equal to ^ the thickness of the channel, which is unequivocally erroneous. 
This is exactly the same method that was used in [21 1 to prove and demonstrate the in- 
stability of a "finger" with a thickness distinct from i with respect to introducing new 
remote poles. The instability of a "finger" with a thickness equal to ^ may be proved 
and demonstrated! This objection was repeatedly made toMark Mineev-Weinstein be- 
fore the publication of his paper [2 IJ . However, no answers can be found there. More- 
over, in our teamwork ll23l . it is shown that for a finite number of poles, any thickness 
"finger" is possible as an asymptotic solution. 

This does not mean, however, that the privileged role of a "finger" no surface ten- 
sion; it only means that the proof is not given in ||2T1 . Let us try to give the correct 
arguments here. The general pole solution (5) in work fSTI is characterised by the real 
parameter a being the sum of the complex parameters ai for poles. The thickness of 
the asymptotic finger is a simple function of a: (Thickness = 1 — ^). The value (a — 1) 
corresponds to a thickness of i. As far as possible, the thickness of the "finger" is be- 
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tween and 1, and the possible a value is in an interval between and 2: (0 < a < 2). 
The value a = 1 corresponding to the finger width ^ is exactly in the middle of this 
interval. What happens to the quite possible initial pole conditions with a outside of 
the limits from to 2? They are "not allowed" because of the already identified finite 
time singularities [23 1. Also, a part of the solutions inside the interval < a < 2 
results in similar finite time singularities. 

Finding exact sufficient and necessary conditions when defining the initial pole 
condition as "not allowed", i.e., singular, is still an open problem. How are these "not 
allowed" initial pole conditions (to be precise, their percentage from the full number 
of possible initial pole conditions corresponding to the given real value a) distributed 
inside of the interval < a < 2? 

From reasons of continuity and symmetry with respect to a = 1 (Fig. [T]!, it is 
possible to conclude that this distribution has a minimum at point a — 1 (thickness ^ !), 
the value that is the most remote from both borders of the interval < a < 2, and that 
the distribution increases to the borders a — 2 and 0, reaching 100 percent for all pole 
solutions outside of these borders, i.e., the thickness ^ is the most probable because 
for this thickness value, the minimal percent of initial conditions potentially capable of 
producing such a thickness value is "not allowed", i.e., results in singularities. A source 
of new poles results in drift of the finger thickness, but this thickness drift is close to 
the most probable and average size equal to i. A similar result is obtained in the case 
of a Saffman-Taylor "finger" with vanishingly small surface tension and with some 
external noise, which was one of the goals of the paper It should be mentioned that 
these formulated arguments are only qualitative and that a strict proof is also necessary. 

3 Asymptotic single Saffman-Taylor "finger" formation 
without surface tension 

In the absence of surface tension, the effect of which is to stabilise the short-wavelength 
perturbations of the interface, the problem of 2D Laplacian growth is described as 
follows: 

{dl + d^y)u = 0. (1) 
M lr(t)= ,9„-u |s= 1 • (2) 

Vn = dnU |r(t) . (3) 

Here, u{x, y; t) is the scalar field mentioned, V{t) is the moving interface, E is a 
fixed external boundary, 9„ is a component of the gradient normal to the boundary (i.e. 
the normal derivative), and u„ is a normal component of the velocity of the front. 

Now, we introduce physical "no-flux" boundary conditions. This means no flux 
occurs across the lateral boundaries of the channel. This requires that the moving 
interface orthogonally intersects the walls of the channel. However, unlike the case 
of periodic boundary conditions, the end points at the two boundaries of the channel 
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do not necessarily have the same vertical coordinate. Nevertheless, this can also be 
considered as a periodic problem in which the period equals twice the width of the 
channel. However, only half of this periodic strip should be considered as the physical 
channel, whereas the second half is its unphysical mirror image. 

Then, we introduce a time-dependent conformal map / from the lower half of a 
"mathematical" plane, = C, + irj, to the domain of the physical plane, z = x + iy, 

where the Laplace equation[T]is defined as ^ — > z. We also require that f{t, ^) w ^ 
for ^ — > C — ioo. Thus, the function z = f{t, C) describes the moving interface. From 
Eqs. ([B, ©, and (O for function f{t, £_) we obtain the Laplacian Growth Equation: 



g^' = 1 k=C-^0 , k k~^oo= 1 • (4) 

Let us look for a solution of Eq. (|4|i in the following form: 



N 

/(^, t) = X^- iT{t) -lY^m log(e^« - e'«'(*)), (5) 
1=1 



N 

a = ^ a; = 1 - A, (6) 

1=1 

where r(f) is some real function of time, ai is a complex constant, ^; = Ci + 
denotes the position of the pole with the number I, and N is the number of poles. 

For our "no-flux" boundary cotextquotedblndition, we must add the condition that 
for every pole ^; = + i^i with ai exists a pole ^/ = — C/ + irji with a7- 

Therefore, we can conclude from this condition for pairs of poles and eq. (|6]l that 
A is a real constant. 

We will prove below that the necessary condition for no finite time singularities for 
a pole solution is 

- 1 < A < 1 , (7) 
Also, for the function F{i£^, t) = if{^, t), for the "no-flux" boundary condition. 



F{ii,t) ^ F{i^,t) (8) 

We want to prove that the final state will be only one finger if no finite time singu- 
larity appears during poles evolutions. 

3.1 Asymptotic behaviour of the poles in the mathematical plane 



This derivation is similar to 11231 . but we assume "no-flux" boundary conditions here (in 
analogy with |24|). The main purpose of this chapter is to investigate the asymptotic 
behaviour of the poles in the mathematical plane. We want to demonstrate that for time 
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t I— ?• oo, all poles go to the two boundary points for no-flux boundary conditions. The 
equation for the interface is 

N 

1=1 

N 

1 - A, -1< A < 1 . (9) 

1=1 

By substitution of Eq. (|9|l in the Laplacian Growth Equation, 

we can find the equations of pole motion (Fig. |2]l: 

Pi = T{t) + (1 - V's^)log h y^a^:log( a^) = const (11) 



and 

AT N 



T = t- 

2 

fe=i 1=1 



where a; = e'^' and Cq is a constant. 
From eqs. ( fTTT i. we can find 



N 

(1 - A)t - log a; 

N N 

Ofett/ log(l - Ofea/) = const . (13) 

fc=i 1=1 

From eqs. (fT2] l and (fTsT l. we can obtain 

/m(^^ a; log ai) — const (14) 



1=1 

and 



1 + A 1 ^ 
t= (^)T + -i?e(^«ilogaO + Ci/2, (15) 

3f the poles, and a; = e'^'^*). 
T M> cxD if t H> cxD and if no finite 

qs. (fTTT i are as follows: 

r + + ^ a; log(l - e'^*'"^''^) = consi, (16) 



where Ci and a; are constants, is the position of the poles, and a; ~ e'^'^* 
In Appendix A, we will prove from eq.(fT2]i that r M> cxd : 

time singularity exists. 

The equations of pole motion that follow from eqs. (fTTT i are as follows 
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a + ^(a;'iog|i-e^(«'-«^) 1 + 

a'l arg(l - e*(«'-«^))) = const, (17) 



arg(l - e'(«'-«^))) = consi, (18) 

where 

6 = 0+*?7/,??/>0. (19) 
a/ = + ia'i' . (20) 

Let us transform 

arg(l -e*(«'"^) = 
arg([l - eKC'-C.)e-(')!+')fc)]) ^ 

arg[l - a,fce''^"-] (21) 

'^/fe=0-a,«a- = e-('"+'"=) (22) 
arg[l — aikd'^'^"'] is a single-valued function of (pik, i.e., 

-|<arg[l-azfee^'^-]<|. (23) 

We multiply eq. ( fTSl l by a'^' and eq. ( fTTI l by a'^, and taking the difference, we 
obtain the following equation: 

a'kCk - a'^T + 
J2{{a'/a',-a';a[)log\l-e^^^'-^^\ + 

l^k 

{a'la'k + a'l'a'l) arg(l - e'(«'-«^))) = const. (24) 

We want to investigate the asymptotic behaviour of poles t !—> oo. 

We have the divergent terms a'^r in this equation. From eq. (l24l). only the term 
log I 1 — e*'^'' I can eliminate this divergence. The necessary condition for this to 
occur is 7]k ^ for t i-^ oo, 1 < k < N. 

We may assume that for t oo, N' groups of poles exist (N' < N) {(pik i— >■ for 
all members of a group). The A^' is currently arbitrary and can even be equal to N. Ni 
is the number of poles in each group, 1 < ^ < N'. 
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For each group, by summation of eqs. (124) over all group poles, we obtain 

«r'cr - «rv + 

5](«"«r - «r"<')iog 1 1 - e'(«r-cr) I + 

l^k 

{af'af + af "af ") arg(l - e'(«f-C))) = const, (25) 



where 

N, 

ar = T.^L (26) 

k 

We have no merging between defined groups for large r, so we investigate the 
motion of poles with this assumption: 

Cr - I > Vr +Vk,l<l,k<N. (28) 
For I ^ k, rif,^ t-^ 0, and (p^^ — Cf^ — CjjJ , we obtain 

log I 1 - e'(«f' -«r) |« log I 1 - e'CCf^-CD | = 

log2 + -logsin2^ (29) 



and 



arg 



(l_e^(C--cr))«arg(l-e'«r-cr)) 



^ + - I . (30) 

We choose n in Eq.(l30]l so that Eq.(l23]l is correct. Substituting these results into eqs. 
(|25] |. we obtain 

Cfe = a^-'CZ - o^V'r + 
E[«"ar-arOloghin^| 

. / qrl qrf , qrfl Qr//\T'lk^ 

+K "fe +af "fc )— ]• (31) 

3.2 Theorem about coalescence of the poles 

From eqs. OTT l. we can conclude the following: 

(i) By summation of eqs. OTI) (or exactly from eq. (fl4l i). we obtain 

ot'k'Ci^ = const . (32) 
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(ii) For I (fii^ |h> 0, 27r, we obtain log | sin ^ cxd, meaning that the poles can 
not pass each other; 

(ill) From (ii), we conclude that <| (p^l \< 27r; 

(iv) From (i) and (iii), (^f,^ ^ cxa is impossible; 

(v) In eq.dSTli. we must compensate for the second divergent term. From (iv) and 
(iii), we can do this only if af^" = a'l — for all /. 



Therefore, from eq. (ISTI . we obtain 



Ni 



^f,:=0, (34) 
7^ , (35) 

Cf = . (36) 

For the asymptotic motion of poles in group N„i, we obtain the following from eqs. 
( [33T l. (O, (|33J, and ([36]l, taking the leading terms in eqs. ([T5j and (fT6] i: 

^ -t, (37) 



A + 1 ' 

= T+> a/ ■ — — — - . (38) 

^ ?7fe +??/ +HCfe -CO 

The solution to these equations is 

r,,=r;^e"4r^*, (39) 



m = ^ife"e -r^'^*, (40) 



a = . (41) 

Therefore, we may conclude that to eliminate the divergent term, we need 

Ni 

ar^Y.<-^^ (42) 
af'(l + A)>0 (43) 

for all /. 
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3.3 The final result 

With the no-flux boundary condition, we have a pair of poles whose condition in eq. 
(l42l l is correct, so all these pairs must merge. Because of the symmetry of the problem, 
these poles can merge only on the boundaries of the channel C, = 0, ±7r. Therefore, 
we obtain two groups of the poles on boundaries. N' = 2, rn = 1, 2, A^i + A^2 ~ N, 
and af + af^' = 1 — A. (In principle, it is possible for some degenerate case of a^^' 
values that eq. (l42T i would be correct for some different groups of poles. However, this 
is a very improbable, rare case.) 

Consequently, we obtain the solution (on two boundaries of the channel Fig. [3]i: 

(1) (i),o ^-^^JTilx* 
Vk - Vk e , (44) 



(1) (i),o ^t^tta* 



Ci'^ = ; (46) 

(2) (2),0 ^^^TTX* 

% - % e °2 , (47) 

^(f ^^,,(2).0e"°r-+^ , (48) 

Cf ^ = ±^ ; (49) 

af'(l + A)>0, (50) 

af'(l + A)>0. (51) 
By summation of eqs. ( |50l ) and ( |50l ) and using eq. |6l we obtain 

(1-A)(1 + A) = 1-A2>0. (52) 



This immediately gives us the formerly formulated condition (|7|i for A 
^ii = ^ — § has an explicit ph 
the moving liquid. We see the 
obtain one finger with width . 



^ii = ^ — § has an explicit physical sense. It is the portion of the channel occupied 
by the moving liquid. We see that for no finite time singularity and for < H> oo, we 

A+l 
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4 Conclusions 

The analytical pole solution for Laplacian growth sometimes yields finite time singu- 
larities. However, an elegant solution of this problem exists. First, we introduce a small 
amount of noise to system. This noise can be considered as a pole flux from infinity. 
Second, for regularisation of the problem, we throw out all new poles that can give a 
finite time singularity. It can be strictly proved that the asymptotic solution for such a 
system is a single finger. Moreover, the qualitative consideration demonstrates that the 
finger equal to ^ of the channel width is statistically stable. Therefore, all properties 
of such a solution are exactly the same as those of the solution with a nonzero surface 
tension under numerical noise. 

Surprisingly, the flame front propagation problem (in spite of exhibiting absolutely 
different physics and mathematical equations for motion) also has analytical pole solu- 
tions and demonstrates the same qualitative behaviour as these solutions in"2ljT6ll . 
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5 Appendix A 



We need to prove that t h- > oo if t h- > c« and if no finite time singularity exists. The 
formula for r is as follows: 

^ N N 

T = t + [--^ ^S^a/ log(l -o^Ta;)] + Co , (53) 

k=l 1=1 

where | a; |< 1 for all 

Let us prove that the second term in this formula is greater than zero: 



^ N N 

■^^^akailog{l -a^ai) = 



k=l 1=1 



^l^l^w^i — 

k=l 1=1 n=l 
oo , N N 



2 ^ n 

n=l k=l 1=1 



oo Af N 

2 E -(E«'(«0")(E > (54) 

n=l 1=1 1=1 



Therefore, the second term in eq. (153) always greater than zero, and consequently, 
T M> oo if i H> cx) for no finite time singularity. 
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Figure Legends 

Fig. 1: The width of the finger is equal to 1 — ^ (The channel width is assumed 
to be equal to 1). The graph for the current a gives the percent of all possible 
solutions resulting in a finite time singularity. The maximum value is equal to 100 
percent and corresponds to a < or a > 2. The minimum is located at the mid- 
dle point a ~ 1 between a = (finger width of 1) and a — 2 (finger width of 0). 
Therefore, at the minimum, the finger width is ^. 

Fig. 2: Geometrical interpretation of the complex constants of motion a'j, = 

and Pk;k^l,...,N. 

Fig. 3: Three consecutive stages of fingering in the Hele-Shaw cell: initial (left), 
intermediate (center), and asymptotic (right). The physical plane z is shown in 
the upper pictures, while the lower pictures depict a distribution of moving poles 
afe(i) in the unit circle < 1 on the mathematical plane uj. The open circle in- 
dicates the repeller, a; = 0, while the solid circle indicates the attractor, cj = 1, of 
poles whose dynamics is given by ( ITT]|T2t . 



REFERENCES 



18 




Figure 2: 



REFERENCES 




Figure 3: 



